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Abstract 

We have measured the probability, n g ^ C c, of a gluon splitting into a charm- 
quark pair using 1.7 million hadronic Z decays collected by the L3 detector. Two 
independent methods have been applied to events with a three-jet topology. One 
method relies on tagging charmed hadrons by identifying a lepton in the lowest 
energy jet. The other method uses a neural network based on global event shape 
parameters. Combining both methods, we measure n g ^ C c = [2.45 ± 0.29 ± 0.53]%. 
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1 Introduction 



The process of a gluon splitting into heavy quark pairs in hadronic Z decays, as shown in 
Fig. [H provides a good test of perturbative QCD [IHlj. The energy scale of heavy quark 
pair production is in the perturbative QCD region, and the process is infrared safe because of 
the natural cutoff provided by the heavy quark mass. It should be mentioned that only the 
process of open heavy quark production is considered here. Measurements of hidden charm 
(J/0) production via QCD processes in Z decays have also been made 

The average number of heavy-quark pairs produced by gluon splitting per hadronic Z decay 
is defined by 

_ _ iV(Z^qqg, g^QQ) 

n g-QQ N ( Z _^ hadrons) ' 1 ' 

where Q stands for a charm or bottom quark. The most recent theoretical predictions jl] 
to leading order in a s , obtained by resumming large leading and next-to- leading logarithmic 
terms to all orders, give n g ^ c5 ~ 2% and ?T g ^ bb - « 0.2%. There are previous experimental 
results for n g _^ cE [OHO] and n g ^ bb - [HUE] • The uncertainty of these measurements contributes 
to the systematic error of electroweak measurements related to heavy quarks and is the biggest 
single source of systematic error in the measurement of _R b = r b b/r had [13J . Therefore, a 
good knowledge of n g ^ cE and ^ g ^ bb is important for precision tests of the Standard Model 
in the heavy-quark sector. In this letter, we present a measurement of n g ^ cE using a high 
statistics sample of hadronic Z decay events collected by L3 [Hj during the years 1994 and 
1995. To identify events with a gluon jet, we first select three-jet events by applying a jet- 
finding algorithm that optimises the yield of g — > cc events. Subsequently, two methods are 
used to identify charm quarks from gluon splitting. The first method, the lepton analysis, 
searches for a lepton from semi-leptonic charm decays in the lowest energy jet. According to 
the JETSET Parton Shower pH] event generator, the latter has about an 80% probability of 
being a gluon jet for the identified three jet sample. The second method, the hadronic event 
shape analysis, uses a neural network technique, with input nodes consisting of several global 
event shape variables, to distinguish g — > QQ events from backgrounds. 



2 Three-Jet Event Selection 

Hadronic events are selected by criteria similiar to the ones used for the measurement of the 
total hadronic cross section jTHj. The number of selected hadronic events is 1.74 million, with 
an estimated background of 0.15% from other processes. To identify the gluon jet, we require a 
three-jet event topology and assume that the lowest energy jet is a gluon jet. The jets are found 
by the JADE algorithm with a y cu t value of 0.03, which maximises the fraction of g — > cc 
events. This cut assigns 63% of g —>■ cc events to the three-jet topology. The jet energies are 
calculated using the relation: 



Ej = E r 



sin ip jk 



sin 



■0jjfe + smip ki + sin0 i; ,-_ 



(2) 



where E cm is the center of mass energy and ipij is f ne angle in space between jets % and j. 
Planar (three jet) events are defined by the condition that the sum of the angles between the 
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three jets is greater than 358°. The jets are labelled in decreasing order of jet energy. To ensure 
that most charged particles are contained in the acceptance of the Silicon Microvertex Detector 
which is used to improve the event signature, the polar angle of the thrust axis, #t, must satisfy 
the condition: | cos^t] < 0.7. With these additional criteria, 430k hadronic events are selected. 

3 Lepton Analysis 

3.1 g — > cc Event Selection and Background Sources 

This method uses events with an electron or a muon candidate in the lowest energy jet to tag 
a charm quark from gluon splitting. Electron candidates are selected by the criteria: 

• An energy cluster with electromagnetic energy greater than 2 GeV, hadronic energy less 
than 2 GeV and polar angle, 9, such that | cos#| < 0.72. 

• The electromagnetic energy cluster must have between 10 and 40 associated BGO crystals 
and more than 95% of the energy in the 9 central crystals. 

• Matching is required in the azimuthal angle and transverse energy between the energy 
cluster and a track in the central track chamber. 

For the muon candidates it is required that: 

• A track is found in the barrel muon system with hits in 2 or more <fi chamber layers and 
1 or more Z- chamber layers. 

• The distance of closest approach to the fill vertex in the transverse plane is less than 100 
mm, as well as less than three times the estimated error due to multiple scattering in the 
calorimeters. 

More details of the electron and muon selection criteria are given in reference |18j . The electron 
energy is required to be between 2 GeV and 6 GeV, while the measured muon momentum 
must be between 3 GeV and 6 GeV; the lepton selection gives 1000 electron and 1287 muon 
candidates. For background and acceptance studies, a Monte Carlo sample of 5.2 million 
hadronic Z decays is used. The events are generated using the JETSET 7.3 generator [15J, 
passed through the full L3 detector simulation and analysed in the same way as the data. 
The following backgrounds, to the g — > cc signal contribute, in decreasing order of importance, 
to the selected lepton samples: 

- Z — > bb and Z — > cc events with hard gluon radiation. In this case, the lowest energy jet 
contains a lepton coming from a semi-leptonic decay of a primary heavy quark. 

- Hadrons misidentified as lepton candidates. 

- Decay in flight of tx^ or K ± into a muon. 

- Dalitz decays: n°,r] — > e + e~j. 

- Photon conversions into electron pairs. 
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- Z -> qqg, g -> bb. 

Monte Carlo studies show that 62% of the selected events with an electron and 47% of the 
events with a muon contain a directly produced heavy quark, usually a b or b, in the lowest 
energy jet. To reduce this background, we use an event discriminant variable [2D] , D, the 
negative logarithm of the probability that all tracks originate from the primary vertex, which 
has small values when all tracks originate from the primary vertex and large values for events 
containing secondary vertices. Applying the cut D < 1.5 rejects 58% of the b-quark background 
while keeping 85% of the g — > cc events. Further discrimination between gluon splitting events 
and directly produced b or b quarks is provided by the invariant mass of the lowest energy jet, 
Mj 3 ; events are rejected if Mj 3 < 6.5 GeV. 

Since all backgrounds are estimated from Monte Carlo events, we have compared several 
of the crucial Monte Carlo distributions with the data. For example, the muon and electron 
momentum, the event discriminant D, and the effective mass of the lowest energy jet, are shown 
in Fig. |2^-d. Good agreement is seen for all these distributions. 

In order to estimate the background from jet misassignment and to check the heavy flavour 
composition in the data and Monte Carlo, we extract the fraction of Z — ► bb events (R^ et ) in 
the three-jet events from data samples using a double-hemisphere tagging method (20]. The 
result found is Rl jet = 0.2017 ± 0.0027 ± 0.0030, where the first error is statistical and the 
second is systematic, mainly due to uncertainties on the hemisphere correlations and on the 
selection efficiencies of the light and charm quarks. This is significantly different from the 
untuned Monte Carlo value of 0.2146 ± 0.0003. The Z — > bb event fraction in the Monte Carlo 
is therefore corrected by reweighting the events according to the measured R^ et value. The 
uncertaintity in the value of R^ et is taken into account in the systematic error. 

Applying the event discriminant and the third-jet invariant mass cut, the ratio between 
data and Monte Carlo predictions for the fraction of selected leptons in the two highest energy 
jets is found to be consistent with unity: 

fDATA rDATA 

= 0.977 ±0.026, ^—^ = 1.025 ±0.053. 

In this background enriched sample good agreement is found between data and Monte Carlo. 
After all selection criteria, 360 electron and 450 muon candidate events are selected in the 
lowest energy jet. The remaining backgrounds are estimated from Monte Carlo as shown in 
Table 1. 



3.2 Results 

After subtracting all background contributions except that due to g — > bb, 51 ± 24 and 75 ± 21 
g — ■> QQ splitting events are found for the muon and electron channels, respectively, where the 
errors are statistical only. The average number of charm quark pairs from gluon splitting per 
hadronic event is then given by: 



Nsei e h Br{b^X£u) _ 

N had -e c -2- Br{c -> X£u) e c Br{c -> X£v) 



_ _ iv se i t ui yu ■—> 



Where e c , e h are the selection efficiencies for g — * cc and g — > bb events, respectively, given in 
Table 2, N se i is the number of events after background subtraction, Nhad is the total number 
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of hadronic Z decays and Br(c — > \£u) and Br(b — > X£z/) are the c and b hadron semileptonic 
branching ratios [21], given in Table 3. Using the efficiencies shown in Table 2, we obtain: 

n g ^ cE = [3.06 ± 1.76 -2.59(n g _ bS - 0.26)]% (muons) 
n g ^ c5 = [3.14 ± 1.14- 3.59(n g ^ bB - 0.26)]% (electrons) 

where the errors are from data statistics only. These results correspond to the weighted average 
of the published values of n g ^ b6 0211111 of [°- 26 ± 0.06]%. 

3.3 Systematic Errors 

The different systematic errors for the lepton analysis are presented in Table 3. There are 
three distinct sources: experimental systematic errors, errors from leptonic branching ratio 
uncertainties and modelling errors that can affect both signal and background. The errors 
labelled 'Monte Carlo Statistics' originate from the limited statistics of the Monte Carlo samples 
of the processes g — > cc and g — > bb used to calculate the selection efficiencies e c and e b . For 
each cut or parameter variation all the efficiencies (including e b ) were recalculated and used in 
Eqn.(3) to obtain the values of 5(rTg_ >C c)/n g _ i . C c presented in Table 3. 

The errors on the lepton detection efficiency are estimated by varying the lepton selection 
cuts by ±10% around their nominal values. The errors on the background simulation are 
estimated using 'lepton enriched' samples of electrons and muons. The systematic errors are 
estimated as the difference between data and Monte Carlo after applying the invariant mass and 
event discriminant cuts to these samples. This difference in the background due to misidentified 
hadrons, photon conversions, Dalitz decays and decays in flight amounts to 2.0% for muons 
and 8.1% for electrons. The 'Track Smearing' error is estimated by smearing tracks to improve 
the data/Monte Carlo comparison for the event discriminant. The error is the difference in 
the gluon splitting rate calculated with, or without, the smeared Monte Carlo. The modelling 
errors are estimated by varying the parameters of the JETSET parton shower model that was 
previously tuned to our data 22 . We vary, within their measured errors, the parameters e& 
and e c in the Peterson fragmentation function [23], the parameter b in the Lund symmetric 
fragmentation function [21] for light quarks, the parameter a q , describing hadron transverse 
momenta and the QCD scale parameter Klla used for the parton shower evolution. Since fully 
simulated events are not available with different values of these parameters, their effects are 
determined from generated events with a detector resolution smeared so as to be consistent 
with the data. The systematic error due to different models of semileptonic decays of charm 
and bottom is also considered. In addition to the JETSET decay model, we estimate the 
uncertainty by also considering the models of Altarelli et aZ.(ACCMM) [22], Isgur et a/.(ISGW) 
and the modified Isgur model (ISGW**) j2S]- The ACCMM model is used for the central values 
of our measurement. 

Combining the muon and electron channels, taking into account the correlated errors arising 
from the leptonic branching ratios and modelling, the overall systematic error is estimated to 
be 21.5%, giving the result: 

n g ^ cE = [3.12 ± 0.96 ± 0.67 - 3.28(n g ^ bb -0.26)]% (leptons) 

where the first error is statistical and the second systematic. 
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4 Event Shape Analysis 



This analysis uses characteristic differences in event shape distributions for the gluon splitting 
to heavy quark process, as compared to the background processes where heavy quarks are 
directly produced by electroweak Z decays. The rather low purity attainable as compared to 
the lepton analysis is largely compensated by data statistics as all hadronic decays of heavy 
quarks are utilised. 

4.1 g — > QQ Event Selection 

In events with a three-jet topology, most of the gluon jets radiated by a primary quark are in 
the two lowest energy jets. Due to the large mass of the heavy quark, a gluon jet containing a 
heavy quark has a larger invariant mass and energy than one containing light quarks. Hence, 
the distribution of the sum of the invariant masses of the two lowest energy jets and the energy 
fraction in a cone around the jet axis are different in a gluon jet with heavy quarks than in 
one with light quarks. Observables, sensitive to the correlation among the three jet momenta, 
can also allow a discrimination between events with a gluon splitting to heavy quark pairs and 
to light quarks or to gluons. In this method, we use three different categories of Monte Carlo 
event samples as follows: 

• 80,000 events containing gluon splitting to charm quark pairs, called the C sample. 

• 8,000 events with gluon splitting to bottom quark pairs, called the B sample. 

• 5.2 million events without the gluon splitting to heavy quark pair process, called the N 
sample. 

A neural network |2Z| has been constructed based on the following five variables: 

• The difference between the sum of the invariant masses of the two lowest energy jets 
and the invariant mass of the highest energy jet, Am = + rri2 — mi, where rrii is the 
effective mass of jet i. 

• The energy in a cone of 8 degree half-angle around the jet axis of the second jet, divided 
by the energy of the jet. 

• Three different Fox- Wolfram moments (n 1 ,n 2 and n 3 ), calculated from the jet momenta, 
sensitive to the global event topology, as described in [28] . 

The neural network has five input nodes, one hidden layer of 10 nodes, and one output node. 
For training, we assume that the samples C and B are the signal events and the N sample is 
background. Subsamples of 10,000 N events, 9,000 C events and 1,000 B events are used for 
training. Since events with a gluon splitting into light quark pairs or gluon pairs are a major 
background source, we use the event discriminant variable, D, to reduce this background. 
Introducing a lower cut on the event discriminant variable in this analysis, the light quark 
background is reduced and the data sample is almost uncorrelated with the one used in the 
lepton analysis. FigO] shows the distribution of the neural network output, O, for g — > QQ and 
for background after applying the cut D > 1.0. Optimising the total error on the signal, we 
choose the cut: O > 0.59. The corresponding selection efficiencies of the data and of the three 
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Monte Carlo samples are listed in Table |U Requiring only the three-jet event selection, the 
purity for the g — > cc events is found to be 2.5%. After the cuts O > 0.59 and D > 1.0, the 
purity is increased to 4.5% and systematic uncertainties are reduced. The purity estimations 
given correspond to the value of 0.15% for n g ^ h ^ used in the JETSET generator. In the region, 
O > 0.59 and 1.0 < D < 1.5, only 35 events from the lepton analysis are found among the 
9,520 selected events. The lepton and event shape analyses can therefore be considered as 
uncorrelated. 

4.2 Results and Systematic Errors 

We extract the probability of a gluon splitting into charm quark pairs using the relation: 

_ _ £ D — £N — {^B — ^v)^ g ^bb , -s 

where Ed, £c, £ b and are the efficiencies for the data, C, B and N samples, respectively. 
The values obtained are listed in Tableland they yield the result: 

n g ^ cE = [2.27 ± 0.30 - 3.86(n g ^ bb - 0.26)]%, 

where the error is from data statistics only, and the dependence on n g ^ bb is shown explicitly. 
This value is in good agreement with that obtained by the lepton analysis. 

The statistical significance of the observed signal is illustrated in Fig 4. The distributions of 
the neural network output for the data and the background Monte Carlo are shown in Fig 4a. 
The relative normalisations are determined according to Eqn.(4) in the region O > 0.59. The 
corresponding g — > QQ signal after background subtraction is shown in Fig 4b as compared to 
the Monte Carlo prediction normalised to the measured value of n g ^ c5 . Although the analysis 
has a low purity, the gluon splitting signal is seen to be large. For the region O > 0.59, it 
amounts to 2064 ± 182 events. The hatched area in Fig 4b shows the estimated g — > bb 
contribution. 

Systematic errors in the efficiencies arise due to imperfect modelling of the event shape 
distributions in the Monte Carlo. We estimate the systematic errors due to several input pa- 
rameters in the JETSET Parton Shower generator by varying their values within the estimated 
errors. For this purpose, fast simulations which take into account the detector resolution have 
been performed for N, C and B samples. Comparing the efficiencies of the new samples, 
generated by variation of one of the parameters, to those obtained with the tuned parameter 
from the reference sample, the changes <5(n g ^ cE )/?T g ^ cE , due to the uncertainties of different 
parameters, are found and listed in Table 

In addition to the tuned fragmentation parameters in JETSET, we consider the following 
sources of systematic error: 

• Monte Carlo statistics. 

• R^ et : We correct the Monte Carlo samples by the measured value of R^ et in three-jet 
events. The changes due to the measurement error on this quantity are assigned as a 
systematic error. 

• Track smearing: This error is estimated in the same way as for the leptonic analysis (see 
section 3.3 above). 
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• Energy calibration: since we reconstruct three-jet events with calorimetric clusters, the 
systematic uncertainty from the energy calibration of the calorimeters is studied. We 
use the Monte Carlo samples of reconstructed three-jet events with energy smearing of 
clusters according to the nominal calorimeter resolution. We repeat the analysis and 
assign the difference as the error. 

Our measurement is also checked by using different values of y cut . With the values of 0.02 and 
0.04 of y cut , we measure the n g ^ cE to be (2.48 ±0.39)% and (2.03 ± 0.27)%, respectively, where 
the error is only statistical. These values are very consistent with that given by the standard 
value y cu t = 0.03. The stability of the result to variation of the event discriminant and neural 
network cuts was also investigated. No systematic deviations were observed. Adding all the 
systematic errors in quadrature, we obtain a total systematic error of 23.7%, yielding the result: 

n g ^ c - c = [2.27 ± 0.30 ± 0.54 - 3.86(rT g _ bb - 0.26)]% (eventshape). 

5 Combined Result and Discussion 

The measurements from the two different methods are now combined in an uncorrelated manner 
giving the result: 

rTg-.cc = [2.45 ±0.29 ±0.53]% , 

where it is assumed that n g _^ bb = [0.26±0.06]%. In the errors quoted, statistical and systematic 
errors (assumed uncorrelated in the two different analysis methods) are combined in quadrature. 
Our result is compared with other measurements and theoretical predictions in Table 6. Good 
agreement is found between the different experimental measurements. The latest resummed 
perturbative QCD calculation [4] agrees, within errors, with our measurement. Our result is 
also in good agreement with the prediction of the ARIADNE and JETSET generators. 
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Category 


Muon channel 


Electron channel 


Observed Events 


450 


360 


Jet misidentification 


153±7 


147±7 


Misidentified hadrons 


188±8 


101±6 


Photon conversions 




2±1 


Dalitz decays 




35±4 


Decays in flight 


58±4 




Total Background 


399±11 


285±10 



Table 1: Summary of the number of selected events and the backgrounds for the lepton analysis. 
The quoted errors are statistical only. 



Category 


Muon channel 


Electron channel 


e c (%) 


0.40 ± 0.04 


0.54 ± 0.05 


e\%) 


0.97 ± 0.10 


1.81 ± 0.14 



Table 2: Summary of selection efficiencies for g — > cc and g — > bb events for the lepton analysis. 
The quoted errors are due to Monte Carlo statistics only. 
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Source 


C / LI \ 1 LL / Crf \ 

5«^cc)/<-cc(%) 


c* / r> \ / o / fry \ 

S(n e R ^)/rf K ^ c - c (%) 


Monte Carlo statistics 


10.1 


9.5 


Lepton efficiency 


6.8 


7.8 


Background simulation 


10.0 


15.0 


Track smearing 


1.9 


3.5 


Total experimental errors 


15.9 


19.7 


Br(b ->^)=0.105± 0.005 


9.8 


4.1 


Br(b -> c -> f)=0.093 ± 0.007 


8.9 


2.1 


Br(c -> f)=0.098± 0.005 


13.7 


9.6 


i? b = 0.2017 ± 0.0040 


4.0 


2.9 


Branching ratio errors 


19.5 


11.0 


A LLA = 0.30 ± 0.03 


4.3 


4.9 


Symmetric parameter o = 0.76±0.08 


2.0 


2.6 


<jq — u.oy±u.uo 


O P. 


O P. 


e b = -0.0056 ±0.0024 


1.6 


1.6 


e c = -0.05 ±0.02 


1.1 


1.1 


Semi-leptonic Decay model 


8.2 


9.3 


Modelling errors 


10.0 


11.3 


Total 


27.1 [15.8] 


25.2 [15.8] 



Table 3: Summary of systematic errors for the muon and electron analysis. The correlated 
error due to branching ratios and QCD modelling is given in square brackets in the last line. 



Event Category 


Efficiency(%) 


Data 
N Sample 
C Sample 
B Sample 


e D = 1.509 ± 0.009 
e N = 1.411 ± 0.006 
e c = 4.401 ± 0.082 
e B = 12.967 ± 0.470 



Table 4: Summary of the selection efficiencies in data and Monte Carlo for the event shape 
analysis. The quoted errors are due to Monte Carlo statistics only. 
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Source 


5(rT g ^cc)/%^cc(%) 


MC statistics 


9.2 


R h = 0.2017 ±0.0040 


14.4 


Track smearing 


2.7 


Energy calibration 


12.5 


A LLA = 0.30 ± 0.03 


4.8 


a q = 0.39 ± 0.03 


7.6 


Symmetric parameter b = 0.76 ± 0.08 


1.7 


e b = -0.0056 ±0.0024 


2.9 


e c = -0.05 ±0.02 


3.8 


Total 


23.7 



Table 5: Summary of systematic errors for the event shape analysis. 



Experimental measurements of n g _> cg 


Expt. 


Ref. 


n g _> C c (%) 


L3 


This paper 


2.45 ±0.29 ±0.53 


OPAL 


8 


3.20 ±0.21 ±0.38 


ALEPH 


9 


3.23 ±0.48 ±0.53 


Theoretical predictions for n g ^ cE 


Model 


Ref. 


n g ^cc (%) 


LO pQCD 


1 


0.607 


Resummed LO pQCD (A) 


3 


1.35 


Resummed LO pQCD (B) 


4 


2.01 


HERWIG 


3 


0.923 


JETSET 


3 


1.70 


ARIADNE 


3 


2.18 



Table 6: Summary of published experimental measurements on n g ^ cg and theoretical pre- 
dictions from perturbative QCD calculations (pQCD) and Monte Carlo models (HERWIG, 
JETSET and ARIADNE). 
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Figure 1: The lowest-order Feynman diagram for a gluon splitting into heavy-quark pairs, 
where Q represents a charm or bottom quark. 
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Figure 2: Distributions of: a) the momentum of the muon candidates and b) the energy of 
electron candidates in the lowest energy jet, c) the event discriminant, D, and d) the invariant 
mass of the lowest energy jet. Data are points with error bars and the open histogram is the 
JETSET Monte Carlo prediction. The shaded histogram shows the JETSET prediction for 

g -»• QQ- 
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Figure 3: Neural network output distribution for g — > QQ (solid line) and background events 
(dashed line). The distributions are normalised to the same area. 
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Figure 4: a) Distribution of the neural network output for data (points), background (his- 
togram) and g — > QQ events (hatched area), b) background-subtracted neural network output 
distribution for data (points) and Monte Carlo g — > QQ events (histogram). The hatched area 
shows the g — > bb contribution. The contributions of g — > cc Monte Carlo events in a) and 
b) are normalised to the measured value of n g ^ cE and the contributions of g — > bb events are 
estimated using the measured value of TT bb- 
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